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(54) LIABLE CARDIAC MUSCLE DIAGNOSING APPARATUS AND LIABLE CARDIAC MUSCLE 
ANALYZING METHOD BY MAGNETIC FIELD MEASUREMENTS 



(57) A magnetic field distribution measurement de- 
vice (1 ) provides a non-contact magnetic measurement 
on a subject's chest at a plurality of coordinates and 
forms therefrom time-series magnetic field distribution 
data. A first arithmetic device (2) in response generates 
image data representing a three-dimensional, intramy- 
ocardial current density distribution. A second arithmetic 
device (3) receives a plurality of tomographic image da- 
ta separately obtained by a tomographic diagnosis ap- 



paratus and processes the data to generate three-di- 
mensional, anatomical image data. A display device (4) 
receives these data and displays on an anatomical im- 
age an image representing an intramyocardial current 
density. Thus in a ventricle a lesioned or viable part of 
myocardium can be readily identified in location, size, 
geometry and level. Furthermore, the anatomical image 
may be replaced with an image representing a normal 
stimulation propagation circuit and serving as a tem- 
plate. 
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Description 

Technical Field 

[0001] The present invention relates generally to via- 
ble myocardium diagnosis apparatuses and analysis 
methods employing a magnetic measurement and par- 
ticularly to those employing a non-contact magnetic 
measurement to non-invasively identify a lesioned part 
of myocardium that might cause an oscillation loop (a 
reentry circuit) of an abnormal electric current in right 
and left ventricles or a viable part of myocardium. 

Background Art 

[0002] Conventionally, electrocardiograms has been 
generally adopted as a technique to diagnose heart dis- 
eases. 

[0003] However, conventional electrocardiography is 
insufficient for example to determine the location, size 
and geometry of a lesioned or viable part of myocardium 
prior to heart surgeries. 

[0004] This is attributed to the fact that electrocardi- 
ography is an indirect measurement methodology. Dif- 
ferent subjects have different tissues existing between 
their hearts and body surfaces, different positional rela- 
tionships between their hearts and other organs and 
bones, their respective hearts having different sizes, a 
different electric conductance for each tissue of their 
bodies, and the like. As such, it has been significantly 
difficult to accurately determine an affected part from in- 
formation obtained from electrocardiography. 
[0005] If a part of myocardium that is lesioned from 
ischemia and a part of myocardium that is physiologi- 
cally viable (hereinafter referred to as "viable myocar- 
dium") coexist in a ventricle, an electrically abnormal 
reentry circuit may be formed in myocardium and for ex- 
ample cause ventricular tachycardia. Accordingly there 
exist a strong demand for accurately, three-dimension- 
ally identifying such a lesioned or viable part of myocar- 
dium. 

[0006] Electrocardiography is useful in diagnosing 
whether myocardium has a part lesioned from ischemia. 
However, it would hardly be so in recognizing the loca- 
tion, size and geometry of a lesioned part of myocardium 
in a heart. 

[0007] In recent years, magnetic resonance imaging 
(MRI), x-ray CT and other medical, tomography appa- 
ratuses have prevailed rapidly. They are useful in exam- 
ining an anatomical structure of the cardiovascular sys- 
tem with high precision. However, they can hardly be 
used to recognize the location, size and geometry of le- 
sioned and viable parts of myocardium. 
[0008] Accordingly, medical institutions currently 
evaluate viable myocardium using myocardial SPECT, 
PET or any other similar tests employing a radioactive 
isotope. Furthermore, evaluating a property of a myo- 
cardial tissue entails an invasive test methodology di- 



rectly obtaining a tissue for example using a catheter or 
performing a surgical operation to sample a myocardial 
tissue (biopsy of myocardium). 
[0009] Of these techniques, myocardial SPECT, PET 

5 or other similar test methodologies using a radioactive 
isotope is used as a method for accurately evaluating a 
survival condition of myocardium before and after angi- 
na pectoris following myocardial infarction, or a coro- 
nary bypass graft operation. These test methodologies, 

10 however, use a radioactive isotope. Accordingly, em- 
ploying them in a controlled area is obligatory and also 
expensive. 

[0010] Furthermore, making a decision after a heart 
transplant as to whether a treatment is necessary in ac- 

15 cordance with the level of a rejection also requires eval- 
uating rejected myocardium, and in addition requires 
conducting a test at relatively small intervals chronolog- 
ically a number of times to see whether a rejection is 
present or absent. Test methodologies employing elec- 

20 trocardiography, echocardiography or the like can be 
employed to capture a variation of a myocardial rejec- 
tion. However, they do not have a precision necessary 
for diagnosis, i.e., as high a precision as necessary for 
determining the level of a rejection, Accordingly after 

25 heart transplants a catheter is used to regularly conduct 
myocardial biopsy, which is a significant invasive on pa- 
tients. Accordingly there has been a strong demand for 
an apparatus capable of employing a non-invasive 
measurement means to three-dimensionally diagnose 

30 the location, size and level of a rejected myocardium af- 
ter heart transplants. 

[0011] In a variety of fields a superconducting quan- 
tum interference device (SQUID) magnetometer has 
been applied. It uses an SQUID capable of detecting 

35 with high sensitivity a magnetic flux of one billionth of 
geomagnetism. In particular, in the field of somatometry, 
which strongly demands non-invasive measurement, as 
described above, an attempt is being made to use a 
SQUID magnetometer to provide a non-contact mag- 

40 netic measurement of human bodies. 

[0012] In particular, the development of thin-film de- 
vice fabrication technology in recent years has allowed 
the development of a DC-SQUID, and an attempt is be- 
ing made to use a SQUID magnetometer to measure a 

45 magnetocardiogram, a distribution of a magnetic field of 
a heart. 

[0013] However, a magnetocardiogram alone cannot 
directly display the location, size, geometry and level of 
a lesioned or viable part of myocardium in a human body 
50 and hardly allows doctors to correctly understand a rel- 
ative, positional relationship of an affected part in a 
heart. 

[0014] Accordingly it has been proposed to visualize 
an intramyocardial, electric current behavior from a 
55 magnetocardiographically represented cardiomagnetic 
field distribution . One such approach adopted is to use 
one or more current dipoles to mimic the source of a 
magnetic field for visualization. However, the entirety of 



2 



3 



EP1 302 160 A1 



4 



myocardium is electrically active with a spread. As such, 
substituting the entirety of myocardium with a single 
electric-current vector cannot provide current density in- 
formation indicating myocardium's electrophysiological 
activity and accordingly it cannot be used to identify the 
size and geometry of a lesioned or viable part of myo- 
cardium. 

[0015] The present invention therefore contemplates 
a viable-myocardium diagnosis apparatus and analysis 
method employing a magnetic measurement, capable 
of employing a non-invasive magnetic measurement to 
obtain data representing an intramyocardial current 
density distribution and used to three-dimensionally 
identify a lesioned or viable part of myocardium safely, 
rapidly and with high precision. 

Disclosure of the Invention 

[0016] In accordance with the present invention an 
apparatus employing a magnetic measurement to diag- 
nose viable myocardium includes a magnetic field dis- 
tribution measurement device, a first arithmetic device, 
a second arithmetic device and a display device. The 
magnetic field distribution measurement device per- 
forms a non-contact magnetic measurement on a sub- 
ject's chest at a plurality of coordinates to obtain a plu- 
rality of time-series magnetic data corresponding to the 
plurality of coordinates, and also uses the plurality of 
time-series magnetic data to generate time-series mag- 
netic field distribution data on the chest. The first arith- 
metic device uses the generated time-series magnetic 
field distribution data to generate time-series, intramy- 
ocardial current density distribution data of the subject. 
The second arithmetic device processes separately pro- 
vided, tomographic, thoracic data of the subject to gen- 
erate data representative of an anatomical image. The 
display device displays an image of the intramyocardial 
current density distribution represented by the data gen- 
erated by the first arithmetic device, as superimposed 
on the anatomical image represented by the data gen- 
erated by the second arithmetic device. Thus a lesioned 
part of myocardium presenting an abnormal current 
density distribution or a physiologically viable part of my- 
ocardium can be identified three-dimensionally. 
[0017] Preferably the first arithmetic device operates 
based on an anatomical or functional factor to divide 
right and left ventricles into any plurality of regions and 
generates time-series data of an intramyocardial current 
density distribution for each the region to allow a le- 
sioned part of myocardium in the right and left ventricles 
to be three-dimensionally identified. 
[0018] More preferably the apparatus further includes 
a database containing information used to determine a 
relationship between an intramyocardial current density 
represented by data generated by the first arithmetic de- 
vice and a myocardial lesion current. 
[001 9] In accordance with the present invention in an- 
other aspect an apparatus employing a magnetic meas- 



urement to diagnose viable myocardium includes a 
magnetic field distribution measurement device, an 
arithmetic device and a display device. The magnetic 
field distribution measurement device performs a non- 

5 contact magnetic measurement on a subject's chest at 
a plurality of coordinates to obtain a plurality of time- 
series magnetic data corresponding to the plurality of 
coordinates, and also uses the plurality of time-series 
magnetic data to generate time-series magnetic field 

10 distribution data on the chest. The arithmetic device us- 
es the generated time-series, magnetic field distribution 
data to generate time-series, intramyocardial current 
density distribution data of the subject. The display de- 
vice uses the data generated by the arithmetic device 

15 to superimpose together an image representing a stim- 
ulation propagation path of the subject extending from 
a sinoatrial node to a bundle of His-Purkinje fiber net- 
work and an image representing an intramyocardial cur- 
rent density distribution and displays the images. Thus 

20 a lesioned part of myocardium presenting an abnormal 
current density distribution or a physiologically viable 
part of myocardium can be three-dimensionally identi- 
fied. 

[0020] Preferably the arithmetic device operates 

25 based on an anatomical or functional factor to divide 
right and left ventricles into any plurality of regions and 
generate time-series data of an intramyocardial current 
density distribution for each the region to allow a le- 
sioned part of myocardium in the right and left ventricles 

30 to be three-dimensionally identified. 

[0021] More preferably the apparatus further includes 
a database containing information used to determine a 
relationship between an intramyocardial current density 
represented by data generated by the arithmetic device 

35 and a myocardial lesion current. 

[0022] In accordance with the present invention in an- 
other aspect a method employing a magnetic measure- 
ment to analyze viable myocardium includes the steps 
of: performing a non-contact magnetic measurement on 

40 a subject's chest at a plurality of coordinates to obtain a 
plurality of time-series magnetic data corresponding to 
the plurality of coordinates and used to generate time- 
series, magnetic field distribution data of the chest and 
generating first data corresponding to time-series, in- 

<5 tramyocardial current density distribution data of the 
subject from the generated time-series magnetic field 
distribution data; processing separately fed, tomograph- 
ic, thoracic image data of the subject to generate second 
data representative of an anatomical image; and dis- 

50 playing an image of the intramyocardial current density 
distribution represented by the first data, as superim- 
posed on the anatomical image represented by the sec- 
ond data, to allow three-dimensional identification of a 
lesioned part of myocardium presenting an abnormal 

55 current density distribution or a physiologically viable 
part of myocardium. 

[0023] Preferably the step of generating the first data 
divides right and left ventricles into any plurality of re- 
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gions based on an anatomical or functional factor and 
generate time-series data of an intramyocardial current 
density distribution for each the region to allow a le- 
sioned part of myocardium in the right and left ventricles 
to be three-dimensionally identified. 5 
[0024] More preferably the method further includes 
the step of determining from information indicating a re- 
lationship between an intramyocardial current density 
represented by the first data and a myocardial lesion 
current whether the three-dimensionally identified part 10 
is a lesioned or viable part of myocardium. 
[0025] In accordance with the present invention in still 
another aspect a method employing a magnetic meas- 
urement to analyze viable myocardium includes the 
steps of: performing a non-contact magnetic measure- *5 
ment on a subject's chest at a plurality of coordinates to 
obtain a plurality of time-series magnetic data corre- 
sponding to the plurality of coordinates and used to gen- 
erate time-series magnetic field distribution data of the 
chest and generating time-series, intramyocardial cur- 20 
rent density distribution data of the subject from the gen- 
erated time-series magnetic field distribution data; and 
using the generated data to superimpose together an 
image representing a stimulation propagation path of 
the subject extending from a sinoatrial node to a bundle 25 
of His-Purkinje fiber network and an image representing 
an intramyocardial current density distribution and dis- 
playing the images to allow three-dimensional identifi- 
cation of a lesioned part of myocardium presenting an 
abnormal current density distribution or a physiological- 30 
ly viable part of myocardium. 

[0026] Preferably the step of generating data divides 
right and left ventricles into any plurality of regions 
based on an anatomical or functional factor and gener- 
ate time-series data of an intramyocardial current den- 35 
sity distribution for each the region to allow a lesioned 
part of myocardium in the right and left ventricles to be 
three-dimensionally identified. 

[0027] More preferably the method further includes 
the step of determining from information indicating a re- 40 
lationship between an intramyocardial current density 
represented by the data and a myocardial lesion current 
whether the three-dimensionally identified part is a le- 
sioned or viable part of myocardium. 
[0028] Thus in accordance with the present invention 45 
an image representing an intramyocardial current den- 
sity distribution obtained through a non-invasive mag- 
netic measurement that is superimposed on an anatom- 
ical image obtained by processing tomographic, thorac- 
ic image data of the same subject undergoing tomogra- so 
phy with a separate, medical diagnosis apparatus can 
be displayed to allow doctors to safely, rapidly and with 
high precision, and three-dimensionally identify a le- 
sioned part of myocardium presenting an abnormal cur- 
rent density distribution or a physiologically viable part 55 
of myocardium. 

[0029] Furthermore in accordance with the present in- 
vention an image representing an intramyocardial cur- 



rent density distribution obtained through a non-invasive 
magnetic measurement that is superimposed on an im- 
age representing a stimulation propagation path of the 
same subject extending from an atrioventricular node to 
a bundle of His-a Purkinje fiber network can be dis- 
played to allow doctors to safely, rapidly and with high 
precision, and three-dimensionally identify a lesioned 
part of myocardium presenting an abnormal current 
density distribution or a physiologically viable part of my- 
ocardium. 

Brief Description of the Drawings 

[0030] In the drawings: 

Fig. 1 is a functional block diagram schematically 
showing a configuration of an apparatus employing 
a magnetic measurement to diagnose viable myo- 
cardium in accordance with the present invention in 
a first embodiment; 

Fig. 2 is a block diagram more specifically showing 
the configuration of the Fig. 1 apparatus; 
Fig. 3 is a block diagram showing a specific config- 
uration of the magnetic field distribution measure- 
ment device shown in Fig. 2; 
Fig. 4 shows by way of example an arrangement of 
a plurality of magnetic sensors on a front side of the 
chest of a subject; 

Fig. 5 represents time-series magnetic data ob- 
tained from the plurality of sensors shown in Fig. 4, 
respectively; 

Fig. 6 schematically illustrates a method of calcu- 
lating current density data from time-series magnet- 
ic data; 

Figs. 7A and 7B show an example of a three-dimen- 
sional, anatomical image displayed on a display de- 
vice 4; 

Fig. 8 tomographically shows a cross section of the 
three-dimensional, anatomical image shown in 
Figs. 7A and 7B; 

Fig. 9 is a flow chart for illustrating an operation of 
the apparatus employing a magnetic measurement 
to diagnose viable myocardium in the first embodi- 
ment; 

Fig. 10 is a flow chart for illustrating an operation of 
the apparatus employing a magnetic measurement 
to diagnose viable myocardium in the first embodi- 
ment in an exemplary variation; 
Fig. 11 schematically shows a manner of dividing a 
ventricle into regions; 

Fig. 12 is a flow chart for illustrating an operation of 
the present apparatus employing a magnetic meas- 
urement to diagnose viable myocardium in the first 
embodiment in another exemplary variation; 
Fig. 13 is a functional block diagram schematically 
showing a configuration of an apparatus employing 
a magnetic measurement to diagnose viable myo- 
cardium in accordance with the present invention in 
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a second embodiment; 
Fig. 14 is a block diagram more specifically showing 
the configuration of the apparatus employing a 
magnetic measurement to diagnose viable myocar- 
dium in the second embodiment shown in Fig. 13; 5 
Figs. 15A and 15B schematically show a normal 
stimulation propagation path in a heart and an elec- 
trocardiograph ically represented waveform; 
Fig. 16 shows an image of a normal stimulation 
propagation path and an abnormal current density *o 
distribution, as displayed by display device 7; 
Figs. 17 and 18 are a flow chart for illustrating first 
and second halves, respectively, of an operation of 
the apparatus employing a magnetic measurement 
to diagnose viable myocardium in the second em- *5 
bodiment; 

Fig. 19 is a flow chart for illustrating a second half 
of an operation of the apparatus employing a mag- 
netic measurement to diagnose viable myocardium 
in the second embodiment in an exemplary varia- 20 
tion; and 

Fig. 20 is a flow chart for illustrating a second half 
of an operation of the apparatus employing a mag- 
netic measurement to diagnose viable myocardium 
in the second embodiment in another exemplary 25 
variation. 

Best Mode for Carrying Out the Invention 

[0031] Hereinafter, the present invention in embodi- 30 
ments will specifically be described with reference to the 
drawings. Note that in the figures, like components are 
denoted by like reference characters and their descrip- 
tions will not be repeated. 

35 

First Embodiment 

[0032] Fig. 1 is a functional block diagram schemati- 
cally showing a configuration of an apparatus in accord- 
ance the present invention in a first embodiment em- 40 
ploying a magnetic measurement to diagnose viable 
myocardium. 

[0033] As shown in Fig. 1 , a magnetic field distribution 
measurement device 1 for example uses a measure- 
ment means such as a SQUID magnetometer, as will be 45 
described hereinafter more specifically, to provide a 
non-contact magnetic measurement on a subject's 
chest at a plurality of coordinates to obtain a plurality of 
time-series magnetic data corresponding to the plurality 
of coordinates. The plurality of time-series magnetic da- 50 
ta are then used to generate and output time-series 
magnetic field distribution data of a magnetic field exist- 
ing on the subject's chest, i.e., of the subject's heart. 
[0034] The cardiac, time-series magnetic field distri- 
bution data provided by magnetic field distribution 55 
measurement device 1 is used by a first arithmetic de- 
vice 2 for example employing a variety of known calcu- 
lation methodologies, described hereinafter, to generate 



and output first data corresponding to time-series, in- 
tramyocardial current density distribution data. 
[0035] Furthermore, magnetic resonance imaging 
(MRI), x-ray, computed tomography (CT), echocardiog- 
raphy, myocardial single photon emission computed to- 
mography (SPECT) or any other similar tomographic di- 
agnosis apparatus is used to separately obtain tomo- 
graphic, thoracic image data (including data of a plurality 
of tomographic images) of the same subject. The data 
are fed to a second arithmetic device 3 and processed 
thereby to generate and output second data represent- 
ing a three-dimensional, anatomical image. 
[0036] The first data is represented in an image. By 
noting a tone of an image representing a current density 
distribution, a lesioned part of myocardium that indi- 
cates an abnormal current density or a viable part of my- 
ocardium can be three-dimensionally identified. 
[0037] Display device 4 superimposes an image rep- 
resenting the intramyocardial current density distribu- 
tion represented by the first data generated by the first 
arithmetic device, on a three-dimensional, anatomical 
image of a subject's chest that is represented by the sec- 
ond data generated by the second arithmetic device 3, 
and displays the same. As a result on an anatomical im- 
age a lesioned or viable part of myocardium can be iden- 
tified three-dimensionally. 

[0038] Note that display device 4 displays only the ex- 
act intramyocardial current density distribution. Accord- 
ingly, a database 5 containing information used to de- 
termine a relationship between an intramyocardial cur- 
rent density distribution and a myocardial lesion current 
is provided and a current density distribution calculated 
by the first arithmetic device 2 is compared with the in- 
formation in database 5 and a result thereof is displayed 
on display device 4 to allow a decision to be more ac- 
curately made as to whether the part having the current 
density of interest is lesioned myocardium (for example 
an ischemic part of myocardium, a rejected myocardium 
and the like) or viable myocardium. 
[0039] Fig. 2 is a block diagram more specifically 
showing the configuration of the apparatus employing a 
magnetic measurement to diagnose viable myocardium 
in accordance with the first embodiment shown in Fig. 1 . 
[0040] As shown in Fig. 2, magnetic field distribution 
measurement device 1 includes in a magnetic shield 
room (MSR) 11 a Dewar structure 13 incorporating a 
SQUID magnetometer and arranged on the chest of a 
subject 12 to provide a non-contact magnetic measure- 
ment, and a magnetic field distribution data operation 
unit 14. 

[0041] In Dewar structure 13 is provided a low-tem- 
perature environment filled with liquid helium to provide 
superconductance, and in the environment is accom- 
modated a SQUID magnetometer configured of a de- 
tector coil formed of superconductor. 
[0042] Fig. 3 is a block diagram more specifically 
showing a SQUID magnetometer 1 5 arranged in an ultra 
low temperature system provided in Dewar structure 13 
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arranged in MSR 11 shown in Fig. 2, and operation unit 
14 arranged in MSR 11 of a normal temperature system. 
[0043] Note that the configuration shown in Fig. 3 is 
that for a single channel for measuring magnetic data of 
a single point on a subject's chest. As will be described 
hereinafter, in the present invention, on a subject's chest 
at a plurality of coordinates a magnetic field is meas- 
ured, i.e., a multi-point, simultaneous magnetic meas- 
urement is provided. Accordingly, MSR 11 of Fig. 2 
would have therein the 1 -channel configuration of Fig. 
3 for each of channels required for a measurement. 
[0044] With reference to Fig. 3, for a single channel a 
SQUID magnetometer generates magnetic data, as de- 
scribed hereinafter. 

[0045] SQUID magnetometer 15 includes a pickup 
coil 16 formed of superconductor for detecting a mag- 
netic field generated from a surface of the chest of a 
subject. When pickup coil 16 captures a magnetic field, 
a current flows and drawn in by a coil 17 to create a 
magnetic field in an Nb shield 20. 
[0046] Consequently, a magnetic field varying linearly 
relative to that created in Nb shield 20 is formed in a 
superconducting loop 18. Voltages of opposite ends of 
superconducting loop 18 are detected by an amplifier of 
operation unit 14 provided in MSR 11 of the normal tem- 
perature system. Operation unit 14 adjusts a current 
flowing through a modulation coil 19 provided in Nb 
shield 20 so that a detected voltage can thus be free of 
variation. 

[0047] More specifically, the detection of an electric 
field of a human body by a SQUID is not a direct meas- 
urement of a magnetic field generated. Rather, a so- 
called a zeropotential method is used to provide a feed- 
back to allow a magnetic field in superconducting loop 
18 to have a constant value (more specifically, a current 
flowing through modulation coil 19 is adjusted to control 
a magnetic field generated in modulation coil 19 so that 
superconducting loop 18 internally, constantly has a 
constant magnetic field) to allow operation unit 14 to 
convert to an electrical signal a magnetic field detected 
at pickup coil 16 and output the signal. Such a feedback 
technique is typically a well known technique referred to 
as a flux locked loop (FLL). 

[0048] Such a SQUID magnetometer 1 5 and its oper- 
ation unit 14 are well known and will not further be de- 
scribed. 

[0049] As has been described previously, the config- 
uration shown in Fig. 3 is that necessary for measuring 
magnetic data for a single channel and outputs an elec- 
trical signal corresponding to time-series magnetic data 
of a magnetic field measured on a front side of the chest 
of a subject at a single point. 

[0050] In the present invention, as has been de- 
scribed previously, a large number of sensors (SQUID 
magnetometers) are arranged on a front side of the 
chest of a subject to measure a magnetic field on the 
front side at multiple points. A magnetic field varies with 
time and for example even during a period correspond- 



ing to a single heart beat a magnetic field that is meas- 
ured at different sites exhibits different variations de- 
pending on the sites. 

[0051] Fig. 4 exemplarily shows an arrangement of a 

5 plurality of sensors (each corresponding to a SQUID 
magnetometer of a single channel) on a front side of the 
chest of a subject. Furthermore, Fig. 5 represents a 
group of time-series magnetic data representing a var- 
iation of a magnetic field for the period of a single heart 

10 beat that is obtained from the respective sensors of Fig. 
4, as corresponding to their respective positions. 
[0052] Magnetic field distribution measurement de- 
vice 1 shown in Fig. 2 outputs a group of time-series 
magnetic data corresponding to a plurality of positions 

15 (coordinates) for measurement, as shown in Fig. 5. If a 
group of time-series magnetic data is captured with a 
particular time noted, it is difficult to graphically (dia- 
grammatically) represent ridges and troughs represent- 
ing a distribution in intensity of a magnetic field present 

20 at a specific time on a front side of a chest to be meas- 
ured, and magnetic field distribution data represented 
in a contour map such as an atmospheric pressure rep- 
resented in a weather chart is accordingly obtained. In 
this sense also, data output from magnetic field distri- 

25 bution measurement device 1 can be captured as time- 
series data of a distribution of a magnetic field on a front 
side of a chest. 

[0053] A group of time-series, magnetic data such as 
output from magnetic field distribution measurement de- 

30 vice 1, i.e., time-series magnetic field distribution data 
are fed to the first arithmetic device 2 of Fig. 2. The first 
arithmetic device 2 functions to obtain from magnetic 
field distribution data of a specific time a density of an 
electric current in the chest that flows at that instant. 

35 [0054] From the time-series magnetic field distribu- 
tion data generated by magnetic field distribution meas- 
urement device 1 the first arithmetic device 2 obtains a 
distribution of a density of a current flowing in a human 
body through a part to be measured (a heart in the 

40 present invention), as described hereinafter. 

[0055] Fig. 6 schematically illustrates a method of ob- 
taining a current density. The method described herein- 
after is intended to indirectly calculate a current that 
should flow through a virtual current sensor arranged at 

45 a specific single part of a human body to be analyzed. 
All of sensors (SQUID magnetometers) on the chest of 
a human body provide time-series magnetic data, which 
can be multiplied by a coefficient and then added togeth- 
er to obtain an output of a current of the virtual sensor 

50 of interest. How the coefficient is obtained is a central 
issue in this arithmetic operation. 
[0056] Hereinafter with reference to Fig. 6 a technique 
used to obtain a current density will be described more 
specifically. Initially on a surface of a human body (a 

55 front side of the chest thereof) there are arranged N 
magnetic sensors in total for the sake of illustration. As 
indicated by a circled number of 1 , a human body (a 
chest, a heart in particular) to be analyzed is regarded 



6 



EP1 302 160 A1 



12 



11 

as a collection of voxels, each in the form of a small 
block. There are provided M voxels in total for the sake 
of illustration. As indicated by a circled number of 2, at 
a single coordinate there is arranged a voxel passing a 
distribution current of mutually orthogonal components, 5 
although a component orthogonal to a plane, as indicat- 
ed by a circled number of 3, is often omitted, since in 
magnetocardiography a magnetic sensor is often ar- 
ranged on a chest in a plane. 
[0057] B|(t) represents time-series magnetic data ob- 
tained from each sensor j and Py represents a spatial 
filter factor of a voxel i corresponding to each sensor 
output Bj(t). 

[0058] If it is assumed that in voxel i there exist a vir- 
tual current sensor and Sj(t) represents a virtual sensor 
output corresponding to a current density obtained from 
the virtual current sensor then S|(t) is defined by the fol- 
lowing expression: 

Si(t) = fp r B.(t). 

j=! 

[0059] As such, if a spatial filter factor Py is determined 
a current density in each voxel i can be obtained and a 
three-dimensional current density distribution can be 
obtained for the entirety of a subject to be analyzed. 
[0060] Spatial filter factor py that is highly sensitive on- 
ly to a distribution current of corresponding voxel i can 
be set by synthetic aperture magnetometry (SAM), mul- 
tiple signal classification (MUSIC) or other similar, vari- 
ous techniques. SAM and MUSIC have been studied 
and developed for example in the fields of radar and so- 
nar and are well known techniques. However, they have 
hitherto been unapplied to cardiomagnetism diagnosis. 
[0061 ] A virtual sensor output of each voxel that is cal- 
culated in real time using a spatial filter factor obtained 
by SAM, MUSIC or any other similar technique, is an 
advantageously highly implemented real-time output. 
[0062] SAM and MUSIC are well known techniques 
and the algorithms using these techniques to obtain a 
spatial filter factor are significantly complicated, and 
they will not be described specifically. SAM is specifical- 
ly described by Robinson SE and Vrba J, "Functional 
Neuroimaging by Synthetic Aperture Magnetometry 
(SAM) m \n Proceedings of the 11th International Confer- 
ence on Biomagnetism, "Reent Advances in 
Biomagnetism, "published by Tohoku University Press, 
1999, pp. 302-305. MUSIC is specifically described by 
Hiroshi Hara and Shinya Kurishiro, "Science ofCerebric 
Magnetism-SQUID Measurement and Medical 
Applications, "published by Ohmsha, January 25, 1997, 
pp. 117-119. 

[0063] Thus the first arithmetic device 2 generates 
from magnetic field distribution data generated by mag- 
netic field distribution measurement device 1 time-se- 
ries data representing a three-dimensional current den- 
sity distribution in a heart to be analyzed and feeds the 



time-series data to display device 4 at one input. 
[0064] The second arithmetic device 3 shown in Fig. 
2 receives image data of a plurality of sliced images (for 
example a dozen of such images obtained at a pitch of 
five millimeters) of the chest of the same subject that 
are previously taken using another tomographic analy- 
sis apparatus (not shown) such as MRI, x-ray CT, 
echocardiography or myocardial SPECT with an elec- 
trocardiography synchronization trigger applied. 
[0065] The second arithmetic device 3 processes (in- 
terpolates) the data of the plurality of sliced images and 
subjects the data to three-dimensional, perspective con- 
version from a predetermined point of view to generate 
second data representing an anatomical image. Thus 
forming a three-dimensional, anatomical image from a 
plurality of sliced images is a well-known technique, for 
example as specifically disclosed in Japanese Patent 
Laying-Open No. 11-128224 and International Publica- 
tion WO 98/15226, and will does not be described spe- 
cifically. 

[0066] Thus the second arithmetic device 3 generates 
the second data representing a three-dimensional, an- 
atomical image of the chest of the same subject in a vi- 
cinity of his/her heart and feeds the second data to dis- 
play device 4 at the other input. 
[0067] Display device 4 of Fig. 2 superimposes on a 
three-dimensional, anatomical image of a subject's 
chest based on the second data from the second arith- 
metic device 3 an image based on the first data from the 
first arithmetic device 2 and representing a three-dimen- 
sional, intra myocardial current density distribution. 
[0068] Figs. 7A and 7B each show a manner of dis- 
playing in real time a three-dimensional current density 
distribution superimposed on a three-dimensional, ana- 
tomical image displayed by display device 4. The figures 
show the current density distribution chronologically 
varying as time elapses. 

[0069] Figs. 7A and 7B each show a three-dimension- 
al image obtained by interpolating approximately five to- 
mographic images obtained for example by slicing a 
subject's chest at a pitch of five millimeters. An actually 
displayed image's depth is difficult to represent by draw- 
ing. In each of Figs. 7A and 7B, a drawing configuring 
each image is represented by a plurality of superim- 
posed drawings. As such, it can thus be determined to 
be a perspective, stereoscopic, anatomical image 
formed by combining a plurality of sliced images. 
[0070] Note that Figs. 7A and 7B each show a tomo- 
graphic image with an upper portion corresponding to a 
front side of a human body and a lower portion to a back 
side thereof. Furthermore Figs. 7A and 7B each show 
the tomographic image, as seen upward (at the legs). 
[0071] In each of Figs. 7A and 7B, a set of circles, as 
indicated by A, represents a three-dimensional current 
density distribution superimposed on a three-dimen- 
sional, anatomical image, and each circle's radius rep- 
resents a magnitude of a current density. Alternatively, 
the magnitude of the current density can be represented 
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by a tone of a specific color on a screen. 
[0072] Furthermore, Fig. 8 displays an extracted tom- 
ographic image obtained at a depth of such a perspec- 
tive, stereoscopic, anatomical image as shown in Figs. 
7A and 7B, with a set of circles, similarly as indicated by 5 
A, representing a current density distribution on the to- 
mographic image. 

[0073] The perspective, three-dimensional, anatomi- 
cal image with the three-dimensional, intramyocardial 
current density distribution displayed thereon allows a 10 
doctor to correctly understand a relative, positional re- 
lationship of the intramyocardial current density distri- 
bution on the anatomical image. In particular, if a dis- 
played current density distribution indicates abnormali- 
ty, the doctor can correctly diagnose the location, size *5 
and geometry of a lesioned or viable part of myocar- 
dium. 

[0074] Fig. 9 is a flow chart representing a method ef- 
fected by the magnetocahogaphic diagnosis apparatus 
of the first embodiment to identify an intramyocardial 20 
current density distribution. 

[0075] In Fig. 9 initially at step S1 magnetic field dis- 
tribution measurement device 1 is used to provide non- 
contact magnetic measurement on the chest of a human 
body at a plurality of coordinates, generate a plurality of 25 
time series data, and record the data if necessary. Note 
that the arithmetic operation/operations performed by 
the first arithmetic device 2 using SAM or MUSIC as 
aforementioned is/are performable on time-series data 
supplied in real time. 30 
[0076] Then at step S2 a plurality of MRI images taken 
with an electrocardiography synchronization trigger ap- 
plied are interpolated by the second arithmetic device 3 
(i.e., subjected to three-dimensional, perspective con- 
version from a predetermined point of view) to obtain a 35 
three-dimensional, anatomical image. 
[0077] Then at step S3 an analysis start time point, 
an analysis end time point and an analysis time interval 
are set to be tg, tg and At, respectively. 
[0078] Then at step S4 analysis start time point tg is <o 
substituted into an analysis time point t to start an anal- 
ysis. Then at step S5 until analysis time point t reaches 
analysis end time point tg the following process is effect- 
ed. 

[0079] More specifically at step S6 the first arithmetic 
device 2 uses SAM or MUSIC to process cardiomagnet- 
ic field distribution data corresponding to a designated 
analysis time point t and obtain intramyocardial current 
density distribution data. 

[0080] Then, at step S7 an effective myocardial cur- 50 
rent density is obtained from an effective volume of a 
myocardial tissue and a calculated current density dis- 
tribution. More specifically, the intramyocardial current 
density distribution data calculated at step S6 is correct- 
ed to be a value closer to the effective myocardial cur- 55 
rent distribution. 

[0081] Then at step S8 display device 4 superimpos- 
es the intramyocardial current density on the anatomical 
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image having been subjected to three-dimensional, per- 
spective conversion from the predetermined point of 
view and displays the same. 

[0082] Then at step S9 At is added to analysis time 
point t. 

[0083] Steps S6-S9 are repeated until a decision is 
made at step S5 that analysis time point t has reached 
analysis end time point tg. When it reaches analysis end 
time point tg display device 4 terminates displaying the 
intramyocardial current density distribution superim- 
posed on the anatomical image. 
[0084] Note that, as has been described previously 
with reference to Fig. 1, information of a lesion at the 
part of myocardium of interest can be obtained from the 
information in database 5, as based on an intramyocar- 
dial current density distribution calculated by the first 
arithmetic device 2. 

[0085] Fig. 1 0 is a flow chart illustrating a process pro- 
vided by an exemplary variation of the first embodiment 
and it is identical to the Fig. 9 flow chart, except the fol- 
lowing: 

[0086] More specifically, at step S10, from the infor- 
mation in database 5 regarding a myocardial lesion cur- 
rent that is based on the corrected myocardial current 
density obtained at step S7 the information of the pres- 
ence/absence of a lesion in myocardium at the part of 
interest and the level of the lesion are extracted and si- 
multaneously displayed on display device 4. 
[0087] The first embodiment in a further exemplary 
variation will now be described with reference to Figs. 
11 and 12. In the first embodiment described above, on 
an anatomical image of a heart a current density distri- 
bution is three-dimensionally, visibly displayed. In this 
exemplary variation, right and left ventricles have myo- 
cardium divided into any plurality of regions based on 
an anatomical or functional factor of the heart and for 
each region a current density distribution is calculated. 
[0088] More specifically, if an average of intramyocar- 
dial current densities can be obtained for each region 
and such averages can be compared to each other, 
whether an electric current flow satisfactorily or not can 
be determined for each region, which is a suitable ma- 
terial allowing doctors to determine whether a part of 
myocardium is a lesioned or viable part. 
[0089] In dividing myocardium into regions, an ana- 
tomical characteristic of a heart is noted to divide right 
and left ventricles into any plurality of regions or an an- 
atomical characteristic of a coronary artery is noted to 
divide a left ventricle into any plurality of regions. Alter- 
natively, for example from an interventricular septum ra- 
dially for every six degrees a density distribution is ob- 
tained to create a profile curve and furthermore create 
a functional representation of voltage and time of elec- 
trical conduction. 

[0090] Fig. 11 schematically shows a manner of divid- 
ing myocardium into regions when the position of a cor- 
onary artery is noted as an anatomical characteristic. In 
the Fig. 11 example, myocardium is divided into a region 
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dominated by a front ramus descendens of a left coro- 
nary artery, that dominated by a left ramus circumflexus, 
and that dominated by a right coronary artery. For each 
of these regions a myocardial current density is calcu- 
lated. 

[0091] Fig. 12 is a flow chart of a process provided by 
a further exemplary variation of the first embodiment 
and it is identical to the Fig. 9 flow chart, except for the 
following: 

[0092] More specifically, at step S11, in a three-di- 
mensional, anatomical image obtained at step S2 myo- 
cardium is divided into regions each dominated by a 
branch of a coronary artery, as shown in Fig. 11 . Then 
at step S12 from the myocardial current density calcu- 
lated at step S7 a myocardial current density is obtained 
for each region provided at step S11 . As such, in addi- 
tion to three-dimensionally displaying a current density 
distribution in accordance with the first embodiment, a 
region of right and left ventricles that includes a lesioned 
part of myocardium can be readily, three-dimensionally 
identified. 

[0093] Thus in accordance with the present invention 
in the first embodiment an image representing an in- 
tramyocardial current density distribution obtained from 
a SQUID magnetometer conducting a non-invasive 
magnetic measurement on a subject's chest that is su- 
perimposed on a three-dimensional, anatomical image 
can be displayed to allow doctors to identify safely, rap- 
idly and with high precision, three-dimensionally identify 
the location, size, geometry and lesion level of a le- 
sioned part of myocardium presenting an abnormal cur- 
rent density distribution or a viable part of myocardium 
to alleviate a burden on patients. Furthermore, a test can 
be conducted without using a radioactive isotope and 
hence continuously without an interval of days. 
[0094] Furthermore, referring to an anatomical or 
functional factor to divide right and left ventricles into a 
plurality of regions and calculating a myocardial current 
density for each region allows three-dimensional diag- 
nosis of a lesioned part of myocardium in the right and 
left ventricles and can thus facilitate considering ap- 
proaches for treatment. In particular, if high frequency 
catheter cauterization is used to provide a treatment, it 
can previously, significantly narrow down a region to be 
electrophysiologically tested using a catheter and a test 
conducted while radioscopy is provided can be conduct- 
ed in a significantly reduced period of time. Consequent- 
ly, doctors and radiographers can avoid significantly 
large annual doses of x-ray exposure. 
[0095] Furthermore, the database that is prepared to 
contain information used to determine a relationship be- 
tween a current density and a myocardial lesion current 
allows a calculated current density distribution to be 
used to provide an accurate diagnosis of an electro- 
physiological property of a tissue of myocardium at the 
part of interest (e.g., an evaluation of viable myocar- 
dium, the presence/absence of a rejection after a heart 
transplant and the level of the rejection, and cardiomy- 



opathy). 

Second Embodiment 

5 [0096] In the first embodiment, forming an anatomical 
image entails obtaining a large number of tomographic 
images of a subject and a test employing MRI, x-ray CT 
or the like is accordingly, previously conducted. This re- 
sults in an increased number of tests and an increased 

10 burden on patients and also an obstacle to a treatment 
directly linked to a test. 

[0097] The present invention in a second embodiment 
can provide a viable-myocardium diagnosis apparatus 
and analysis method employing a magnetic measure- 
is ment capable of eliminating the formation of an anatom- 
ical image to conduct a reduced number of tests and 
carry out a diagnosis and a test such that they are di- 
rectly linked. 

[0098] Fig. 1 3 is a functional block diagram schemat- 
20 ically showing a configuration of the viable-myocardium 
diagnosis apparatus employing magnetic measurement 
in the second embodiment. 

[0099] With reference to Fig. 1 3, magnetic field distri- 
bution measurement device 1 will not be described as 

25 it has been described in the first embodiment. 

[01 00] Magnetic field distribution measurement appa- 
ratus 1 generates time-series, magnetic field distribution 
data and outputs the data to an arithmetic device 6 
which in turn uses the received time-series, magnetic 

30 field distribution data and employs the aforementioned 
SAM, MUSIC or any other similar calculation technique 
to generate data representing an intramyocardial, three- 
dimensional current density distribution. Arithmetic de- 
vice 6 then uses the generated three-dimensional cur- 

35 rent density distribution data to superimpose data rep- 
resenting an intracardiac excitation (stimulation) propa- 
gation path of a period corresponding to that from an 
electrocardiographically represented P wave to an elec- 
trocardiographically represented QRS complex and da- 

40 ta representing a current density distribution on each 
other and outputs the same to display device 7. 
[0101] Display device 7 superimposes an image rep- 
resenting the intramyocardial current density distribu- 
tion represented by the data generated by arithmetic de- 

45 vice 6, on a three-dimensional image of the excitation 
propagation path also obtained by arithmetic device 6 
and corresponding to the period from the P wave to the 
QRS complex, and displays the same. Consequently, 
such an anatomical image as used in the first embodi- 

50 ment can be dispensed with to three-dimensionally 
identify an abnormal current density distribution in my- 
ocardium of a ventricle. 

[0102] Note that as well as in the first embodiment, 
there is provided database 5 containing information 
55 used to determine a relationship between an intramyo- 
cardial current density distribution and a myocardial le- 
sion current. 

[0103] Fig. 14 is a block diagram more specifically 
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showing the configuration of the viable-myocardium di- 
agnosis apparatus employing magnetic measurement 
in the second embodiment shown in Fig. 13. 
[0104] With reference to Fig. 14, magnetic field distri- 
bution measurement apparatus 1 will not be described 5 
as it is identical to that described with reference to Figs. 
2 and 3. 

[0105] Magnetic field distribution measurement de- 
vice 1 outputs time-series, magnetic field distribution da- 
ta and outputs the data to arithmetic device 6 shown in 
Fig. 14. Arithmetic device 6 uses SAM, MUSIC or any 
other similar technique, as described with reference to 
Fig. 6, to convert the time-series, magnetic field distri- 
bution data to time-series, current density distribution 
data. 

[0106] Subject 12 has his/her electrocardiogram re- 
corded by an electrocardiograph 21 to allow measured 
electrocardiographic waveform data of subject 12 to be 
fed to arithmetic device 6. 

[0107] Note herein that if the electrocardiographically 
represented waveform and a generated current density 
distribution can be correlated, the electrocardiogram 
and an event occurring in the heart can also be corre- 
lated. 

[01 08] Reference will now be made to Fig. 1 5A sche- 
matically representing a normal stimulation propagation 
path in a heart and Fig. 15B representing an electrocar- 
diographically represented waveform for a single heart 
beat. 

[0109] With reference to Figs. 15A and 15B, a si- 
noatrial node functions as a pacemaker determining a 
heart beat and it fires at predetermined intervals (a tim- 
ing of a P wave of an electrocardiogram) to generate a 
pulse. This pulse is transmitted through a specific stim- 
ulation propagation path to an atrioventricular node and 
therein after a period of time elapses a pulse is trans- 
mitted through a bundle of His-Purkinje fiber network to 
an underlying a ventricle and myocardial contraction 
erupts. This propagation of a stimulation from the bundle 
of His to Purkinje fiber network corresponds to the peri- 
od of the QRS complex in the electrocardiogram (an is- 
ovolumetric contraction time). 

[0110] As such, by analyzing cardiomagnetism relat- 
ed to the period from the P wave to the QRS complex, 
i.e., an intra myocardial current density distribution, 
arithmetic device 5 generates image data representing 
a stimulation propagation path serving as a normal 
route, as shown in Fig. 15A. 

[0111] An image of a stimulation propagation path, 
such as shown in Fig. 15A, can be used in place of the 
anatomical image used in the first embodiment, as a 
template displayed. More specifically, a three-dimen- 
sional, anatomical image such as described in the first 
embodiment can be dispensed with if a stimulation prop- 
agation path of a normal route, as shown in Fig. 15A, is 
displayed, since a lesioned part of myocardium present 
in a vicinity thereof that presents an abnormal current 
density distribution or a viable part of myocardium would 



be readily, anatomically correlated and its location, size 
and geometry would be identified by doctors. 
[0112] Arithmetic device 6 shown in Fig. 14 generates 
data representing a generated current density distribu- 
tion, superimposed on a displaying of a stimulation prop- 
agation circuit as a template, such as described above. 
As has been described previously, noting an image rep- 
resenting a current density distribution allows a lesioned 
or viable part of myocardium in a ventricle to be found, 
and such image data can be combined with the afore- 
mentioned template image data and fed to display de- 
vice 1. 

[0113] Display device 7 shown in Fig. 14 uses the data 
received from arithmetic device 6 to display an image 
representing a current density distribution, as superim- 
posed on a normal stimulation propagation circuit serv- 
ing as a template. Thus doctors can anatomically corre- 
late a lesioned or viable part of myocardium. 
[0114] Fig. 16 exemplarily shows a screen actually 
displayed by display device 7. It displays an image of a 
part presenting an abnormal current density distribution, 
superimposed on a normal stimulation propagation cir- 
cuit serving as a template. 

[0115] A doctor would be able to use a positional re- 
lationship of the current density distribution relative to 
the normal stimulation propagation path serving as a 
template, as shown in Fig. 16, to readily provide an an- 
atomical correlation and identify the location, size, ge- 
ometry and level of a lesioned or viable part of myocar- 
dium in right and left ventricles. 
[01 1 6] Figs. 1 7and 1 8 are a flow chart representing a 
method effected by the apparatus of the second embod- 
iment to electrophysiologically analyze a viable myocar- 
dium. 

[0117] Initially, with reference to Fig. 17, at step S21 
magnetic field distribution measurement device 1 is 
used to provide a non-contact magnetic measurement 
on the chest of a human body at a plurality of coordi- 
nates to generate and record a plurality of time-series 
magnetic data. 

[0118] Then at step S22 an analysis start time point 
is set to correspond to an electrocardiographically rep- 
resented P-wave start time point tgp, an analysis end 
time point to an electrocardiographically represented 
QRS-complex end time point f^Rg, and an analysis 
time interval to At. 

[0119] Then at step S23 time point tgp is substituted 
into analysis time point t. 

[01 20] Then at step S24 until an analysis time reaches 
time point t eQRS the following steps S25-S27 are repeat- 
ed. 

[01 21 ] More specifically at step S25 arithmetic device 
6 uses SAM or MUSIC to process cardiomagnetic field 
distribution data of a designated analysis time point t to 
generate intramyocardial current density distribution da- 
ta. 

[0122] Then at step S26 the intramyocardial current 
density distribution data having been subjected to three- 
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dimensional, perspective conversion from a predeter- 
mined point of view is displayed in an image. 
[0123] Then at step S27 At is added to analysis time 
point t and the process returns to step S24 and a deci- 
sion is made as to whether time point te QRS has been 
reached. If so then it means that there has been ob- 
tained image data representing a stimulation propaga- 
tion path corresponding to a normal route, as shown in 
Fig. 15A, as corresponding in an electrocardiographi- 
cally represented waveform to the period from the P 
wave to the QRS complex. 

[0124] Then the process proceeds with step S28, 
shown in Fig. 18, and an analysis start time point, an 
analysis end time point and an analysis time interval are 
set to be tg, ^ and At, respectively. 
[0125] Then at step S29 analysis start time point tg is 
substituted into analysis time point t. 
[0126] Then at step S30 until a decision is made that 
analysis time point t has reached analysis end time point 
i e the following steps S31-S34 are effected. 
[0127] More specifically, at step S3 1 arithmetic device 
6 uses SAM or MUSIC to process cardiomagnetic field 
distribution data of a designated analysis time point t to 
generate intra myocardial current density distribution da- 
ta. 

[0128] Then at step S32 an effective myocardial cur- 
rent density is obtained from the effective volume of a 
myocardial tissue and a calculated current density dis- 
tribution. More specifically, the current density distribu- 
tion data calculated at step S31 is corrected to have a 
value closer to the effective myocardial current density. 
[0129] Then at step S33 the intramyocardial current 
density data is superimposed on an image of a normal 
stimulation propagation circuit having been subjected to 
three-dimensional, perspective conversion from a pre- 
determined point of view, and it is thus displayed. 
[01 30] Furthermore at step S34 At is added to analy- 
sis time point t and the process returns to step S30 and 
a decision is made as to whether analysis end time point 
t e has been reached. Thus, data representing the in- 
tramyocardial current density distribution is superim- 
posed on an image of a normal stimulation propagation 
path (Fig. 15A) obtained through the Fig. 17 flow chart, 
and it is thus displayed. 

[0131] Note that as well as in the first embodiment, 
from database 5 of Fig. 13 information can be obtained 
regarding a lesion in myocardium at a lesioned part. 
[01 32] Fig. 1 9 is a flow chart of a process of an exem- 
plary variation of the second embodiment and it is iden- 
tical to the Fig. 18 flow chart, except the following: 
[0133] More specifically, at step S35, from the infor- 
mation in database 5 regarding a myocardial lesion cur- 
rent that is based on the corrected myocardial current 
distribution obtained at step S32 the information of the 
presence/absence of a lesion in myocardium at the part 
of interest and the level of the lesion are extracted and 
simultaneously displayed on display device 7. 
[0134] Then, as well as in the first embodiment, in the 



present embodiment right and left ventricles may also 
have their myocardium divided into any plurality of re- 
gions based on an anatomical or functional, factor of the 
heart and for each region a current density distribution 
5 may be calculated. The myocardium is divided into re- 
gions in the manner as described in the first embodi- 
ment. 

[01 35] Fig. 20 is a flow chart of a process provided by 
the variation of the second embodiment, as described 

10 above, and it is identical to the Fig. 18 flow chart. 

[0136] More specifically, at step S36 the position of 
an interventricular septum is clarified from information 
of a stimulation propagation path extending from a bun- 
dle of His to immediately before a Purkinje fiber network 

15 that is obtained by arithmetic device 6. Accordingly, a 
ventricle is divided into a plurality of regions at a prede- 
termined angle with an atrioventricular node serving as 
a reference. Then, at step S37, from the effective vol- 
ume and current density distribution of each region of a 

20 thus divided tissue a myocardial current density calcu- 
lated for each region can be obtained. 
[01 37] Thus in accordance with the present invention 
in the second embodiment an image representing an in- 
tramyocardial current density distribution obtained from 

25 a SQUID magnetometer obtaining a non-invasive mag- 
netic measurement on a subject's chest that is superim- 
posed on a normal stimulation propagation path serving 
as a template can be displayed to be eliminate the ne- 
cessity of superimposing it on a different, anatomical im- 

30 age to allow doctors to safely, rapidly and with high pre- 
cision, three-dimensionally identify the location, size, 
geometry and level of a lesioned part of myocardium ex- 
hibiting an abnormal current density distribution, or a vi- 
able part of myocardium. Electrophysiolographical func- 

35 tion in the myocardium can thus be non-invasively, rap- 
idly and safely tested. An invasion on patients can thus 
significantly be alleviated. Furthermore, a previous test 
conducted to obtain an anatomical image can be elimi- 
nated. 

40 [0138] Furthermore, a test can be conducted without 
using a radioactive isotope. This can reduce the cost 
and the test can be conducted continuously without an 
interval of days. 

[0139] Furthermore, referring to an anatomical or 
45 functional factor to divide right and left ventricles into a 
plurality of regions and calculating a myocardial current 
density for each region allows three-dimensional diag- 
nosis of a lesioned part of myocardium in the right and 
left ventricles and can thus facilitate considering ap- 
50 proaches for treatment. In particular, if high frequency 
catheter cauterization is used to provide a treatment, it 
can previously, significantly narrow down a region to be 
electrophysiologically tested using a catheter and a test 
conducted while radioscopy is provided can be conduct- 
55 ed in a significantly reduced period of time. Consequent- 
ly, doctors and radiographers can avoid significantly 
large annual doses of x-ray exposure. 
[0140] Furthermore, the database that is used to de- 
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termine a relationship between a current density and a 
myocardial lesion current allows an accurate diagnosis 
of an electrophysiological function of myocardium (e.g., 
an evaluation of viable myocardium, the presence/ab- 
sence of a rejection after a heart transplant and the level 
the rejection, and cardiomyopathy). 
[0141] Thus in the present invention an intramyocar- 
dial current density distribution obtained from a non-in- 
vasive magnetic measurement performed on a patient's 
chest can be visibly displayed on a three-dimensional, 
anatomical image to allow three-dimensional identifica- 
tion of the location, size, geometry and level of a le- 
sioned part of myocardium exhibiting an abnormal cur- 
rent density distribution or a viable part of myocardium. 
This allows non-invasive diagnosis of the lesioned or vi- 
able part of myocardium. As such, a test can be con- 
ducted rapidly and safely without imposing a significant 
burden on patients. Furthermore, a test can be conduct- 
ed without using a radioactive isotope. As such, it can 
be conducted continuously, which can effectively reduce 
the cost. 

[0142] Furthermore, referring to an anatomical or 
functional factor to divide right and left ventricles into a 
plurality of regions and calculating a myocardial current 
density for each region allows three-dimensional diag- 
nosis of a lesioned or viable part of myocardium. In par- 
ticular, if high frequency catheter cauterization is used 
to provide a treatment, it can previously, significantly 
narrow down a region to be electrophysiologically tested 
using a catheter. Consequently, doctors and radiogra- 
phers can effectively avoid significantly large doses of 
x-ray exposure. 

[0143] Furthermore, the database that is prepared to 
contain information used to determine a relationship be- 
tween a current density and a myocardial lesion current 
uniquely effectively allows an accurate diagnosis of an 
electrophysiological property of a tissue of myocardium. 
[0144] In accordance with the present invention, in 
still another aspect, superimposing a current density 
distribution on a normal stimulation propagation circuit 
in the same subject extending from a sinoatrial node to 
a bundle of His-a Purkinje fiber network, and three-di- 
mensionally displaying the same can eliminate the ne- 
cessity of obtaining an anatomical image to three-di- 
mensionally identify the location, size, geometry and 
level of a lesioned or viable part of myocardium in a ven- 
tricle. In other words, a test conducted to obtain an an- 
atomical image can be eliminated to effectively provide 
a more cost-effective diagnosis. 

Industrial Applicability 

[0145] Thus the present invention can provide a via- 
ble-myocardium diagnosis apparatus and analysis 
method employing a magnetic measurement, capable 
of three-dimensionally identifying the location, size, ge- 
ometry and level of a lesioned part of myocardium ex- 
hibiting an abnormal current density distribution or a vi- 



able part of myocardium. It is thus useful in non-invasive 
diagnosis of a lesioned or viable part of myocardium and 
furthermore in a treatment employing high frequency 
catheter cauterization. 



Claims 

1 . An apparatus employing a magnetic measurement 
10 to diagnose viable myocardium, comprising: 

a magnetic field distribution measurement de- 
vice (1) performing a non-contact magnetic 
measurement on a subject's chest at a plurality 
15 of coordinates to obtain a plurality of time-se- 

ries magnetic data corresponding to said plu- 
rality of coordinates, and also using said plural- 
ity of time-series magnetic data to generate 
time-series, magnetic field distribution data on 
20 said chest; 

a first arithmetic device (2) using said generat- 
ed time-series, magnetic field distribution data 
to generate time-series, intramyocardial cur- 
rent density distribution data of said subject; 
25 a second arithmetic device (3) processing sep- 

arately provided, tomographic, thoracic data of 
said subject to generate data representative of 
an anatomical image; and 
a display device (4) displaying an image of said 
30 intramyocardial, current density distribution 

represented by said data generated by said first 
arithmetic device, as superimposed on said an- 
atomical image represented by said data gen- 
erated by said second arithmetic device, there- 
35 by capable of three-dimensionally identifying a 

lesioned part of myocardium presenting an ab- 
normal current density distribution or a physio- 
logically viable part of myocardium. 

40 2. The apparatus of claim 1 , wherein said first arith- 
metic device operates based on an anatomical or 
functional factor to divide right and left ventricles in- 
to any plurality of regions and generates time-series 
data of an intramyocardial current density distribu- 
45 tion for each said region to allow a lesioned part of 
myocardium in the right and left ventricles to be 
three-dimensionally identified. 

3. The apparatus of claim 1, further comprising a da- 
50 tabase (5) containing information used to determine 
a relationship between an intramyocardial current 
density represented by data generated by said first 
arithmetic device and a myocardial lesion current. 

55 4. An apparatus employing a magnetic measurement 
to diagnose viable myocardium, comprising: 

a magnetic field distribution measurement de- 
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vice (1) performing a non-contact magnetic 
measurement on a subject's chest at a plurality 
of coordinates to obtain a plurality of time-se- 
ries magnetic data corresponding to said plu- 
rality of coordinates, and also using said plural- 5 
ity of time-series magnetic data to generate 
time-series magnetic field distribution data on 
said chest; 

an arithmetic device (6) using said generated 
time-series, magnetic field distribution data to 10 
generate time-series, intramyocardial current 
density distribution data of said subject; and 
a display device (7) using the data generated 
by said arithmetic device to superimpose to- 
gether an image representing a stimulation 15 
propagation path of said subject extending from 
a sinoatrial node to a bundle of His-Purkinje fib- 
er network and an image representing an in- 
tramyocardial current density distribution and 
display said images, thereby capable of three- 20 
dimensionally identifying a lesioned part of my- 
ocardium presenting an abnormal current den- 
sity distribution or a physiologically viable part 
of myocardium. 

25 

5. The apparatus of claim 4, wherein said arithmetic 
device operates based on an anatomical or func- 
tional factor to divide right and left ventricles into any 
plurality of regions and generate time-series data of 

an intramyocardial current density distribution for 30 
each said region to allow a lesioned part of myocar- 
dium in the right and left ventricles to be three-di- 
mensionally identified. 

6. The apparatus of claim 4, further comprising a da- 35 
tabase (5) containing information used to determine 

a relationship between an intramyocardial current 
density represented by data generated by said 
arithmetic device and a myocardial lesion current. 

40 

7. A method employing a magnetic measurement to 
analyze viable myocardium, comprising the steps 
of: 

performing a non-contact magnetic measure- 45 
ment on a subject's chest at a plurality of coor- 
dinates to obtain a plurality of time-series mag- 
netic data corresponding to said plurality of co- 
ordinates and used to generate time-series, 
magnetic field distribution data of said chest 50 
and generating first data corresponding to time- 
series, intramyocardial current density distribu- 
tion data of said subject from the generated 
time-series magnetic field distribution data; 
processing separately fed, tomographic, tho- 55 
racic image data of said subject to generate 
second data representative of an anatomical 
image; and 



displaying an image of said intramyocardial 
current density distribution represented by said 
first data, as superimposed on said anatomical 
image represented by said second data, to al- 
low three-dimensional identification of a le- 
sioned part of myocardium presenting an ab- 
normal current density distribution or a physio- 
logically viable part of myocardium. 

8. The method of claim 7, wherein the step of gener- 
ating the first data divides right and left ventricles 
into any plurality of regions based on an anatomical 
or functional factor and generates time-series data 
of an intramyocardial current density distribution for 
each said region to allow a lesioned part of myocar- 
dium in the right and left ventricles to be three-di- 
mensionally identified. 

9. The method of claim 7, further comprising the step 
of determining from information indicating a rela- 
tionship between an intramyocardial current density 
represented by said first data and a myocardial le- 
sion current whether said three-dimensionally iden- 
tified part is a lesioned or viable part of myocardium. 

10. A method employing a magnetic measurement to 
analyze viable myocardium, comprising the steps 
of: 

performing a non-contact magnetic measure- 
ment on a subject's chest at a plurality of coor- 
dinates to obtain a plurality of time-series mag- 
netic data corresponding to said plurality of co- 
ordinates and used to generate time-series 
magnetic field distribution data of said chest 
and generating time-series, intramyocardial 
current density distribution data of said subject 
from the generated time-series magnetic field 
distribution data; and 

using said generated data to superimpose to- 
gether an image representing a stimulation 
propagation path of said subject extending from 
a sinoatrial node to a bundle of His-Purkinje fib- 
er network and an image representing an in- 
tramyocardial current density distribution and 
displaying said images to allow three-dimen- 
sional identification of a lesioned part of myo- 
cardium presenting an abnormal current densi- 
ty distribution or a physiologically viable part of 
myocardium. 

1 1 . The method of claim 1 0, wherein the step of gener- 
ating data divides right and left ventricles into any 
plurality of regions based on an anatomical or func- 
tional factor and generates time-series data of an 
intramyocardial current density distribution for each 
said region to allow a lesioned part of myocardium 
in the right and left ventricles to be three-dimension- 
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aily identified. 

12. The method of claim 10, further comprising the step 
of determining from information indicating a rela- 
tionship between an intramyocardial current density 
represented by said data and a myocardial lesion 
current whether said three-dimensionally identified 
part is a lesioned or viable part of myocardium. 
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( START ) 



EFFECT A NON-CONTACT MAGNETIC MEASUREMENT ON THE CHEST 
OF A HUMAN BODY AT A PLURALITY OF COORDINATES & RECORD 
A PLURALITY OF TIME-SERIES DATA 



SI 



OBTAIN A PLURALITY OF MRI IMAGES WITH AN 
ELECTROCARDIOGRAPHY SYNCHRONIZATION TRIGGER APPLIED 
THERETO & INTERPOLATE THE MRI IMAGES TO OBTAIN A 
THREE-DIMENSIONAL, ANATOMICAL IMAGE 



SET AN ANALYSIS START TIME POINT, AN ANALYSIS END TIME 
POINT & AN ANALYSIS TIME INTERVAL TO BE ts. te &At. 
RESPECTIVELY 



SUBSTITUTE ts INTO AN ANALYSIS TIME POINT t 
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S3 



S4 




( END ) 



PROCESS CARDIOMAGNETIC FIELD DISTRIBUTION DATA OF A 
DESIGNATED ANALYSIS TIME POINT t BY SAM (OR MUSIC) TO 
OBTAIN INTRAMYOCARDIAL CURRENT DENSITY DISTRIBUTION 
INFORMATION 
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OBTAIN (CORRECT) A MYOCARDIAL CURRENT DENSITY FROM A 
TISSUE'S EFFECTIVE VOLUME & A CURRENT DENSITY 
DISTRIBUTION 



I 



SUPERIMPOSE THE INTRAMYOCARDIAL CURRENT DENSITY 
DISTRIBUTION INFORMATION ON AN IMAGE HAVING UNDERGONE 
THREE-DIMENSIONAL, PERSPECTIVE CONVERSION FROM A 
PREDETERMINED POINT OF VIEW & DISPLAY THE SAME 
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